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Abstract 

Continuous interpolation of real-valued data is characterized by piece- 
wise monotone functions on a compact metric space. Topological total 
variation of piecewise monotone function / : X —>■ R is a homeomorphism- 
invariant generalization of ID total variation. A varilet basis is an or¬ 
thonormal collection of piecewise monotone functions {gi \ i = !...«}, 
called varilets, such that every linear combination '^Oigi {m G R) has 
topological total variation 57 I®*I- ^ varilet transform for / is a varilet 
basis for which / = 57 Filtered versions of / result from altering the 
coefficients Oi. 
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1 Introduction 


Topology has proved to be a powerful tool for discovery of the structure of data 
mi [5]. Typically, data is viewed as a finite sample of a continuous function. The 
function becomes the object of study, often via a topological representation such 
as the contour tree [S], Reeb graph m, Morse-Smale complex [3, persistence 
diagram mi or persistence barcode [6]. 

This paper presents new results for topological analysis of continuous functions. 
Whereas computational topology often utilizes algebraic topology mi , we ex¬ 
ploit the older field of analytic topology EH- 

Our starting point is real-valued continuous functions on compact metric spaces. 
We introduce the broad category of piecewise monotone functions, well-suited to 
data interpolation. We introduce a new measure on functions - topological total 
variation - and we introduce varilets, an orthonormal basis which additively 
decomposes both the function and its topological total variation. 

The varilet transform maps a function to a varilet basis by means of a lens pa¬ 
rameter. A lens is collection of upper and lower level set components, providing 
a multiresolution view of the function. 

We state a mathematical algorithm - the Varilet Transform Algorithm - and 
prove its output correct. 

We proceed with an overview, followed by discussion of related work and iden¬ 
tification of contributions. Two sections then develop the analytic topology, the 
first defining the Varilet Transform Algorithm, and the second validating its 
correctness. 

1.1 Overview 

Analytic topology m. continuum theory [2zi ng eg E] and dynamic topology 
PSI Eg had their heyday in the mid-twentieth century. The Eilenberg-Whyburn 
monotone-light factorization uni Eg is a powerful result concerning functions 
on compact metric spaces. 

This paper uses the monotone-light factorization as the foundation for topolog¬ 
ical analysis of real-valued functions. Varilets are an elementary application of 
analytic topology. 

Continuous function f : X ^ Y is monotone when f~^y is connected for all 
y GY] ] thus f~^ carries connected sets to connected sets. / is light when f~^y 
is totally disconnected for all y G Y. The monotone-light factorization EH Eg 
states that there exists a unique compact metric space M - called /’s middle 
space - such that / = A o where p. ■. X ^ M is monotone and A : M —>• V is 
light. 

The middle space M of function f : X ^ Y is the quotient of the domain that 
identifies all points which, for some y G Y, lie in the same connected component 
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of / ^y- /’s monotone factor fjL is the quotient map from X to M\ f’s light 
factor A assigns each point p G M the value G Y. Thus f = Xo p. 

Piecewise monotonicity is defined in terms of the middle space. The middle 
space M of piecewise monotone scalar field / : X —>■ IR is a graph having finitely 
many vertices and edges, with each edge having the topology of a closed real 
interval. M is identical to /’s Reeb graph [2T]. The monotone factor p locates 
extrema, saddles and contours in the domain X ; the light factor A provides their 
numerical values. 

The monotone-light factorization enables one to proceed by first defining con¬ 
structions using middle space M and light factor A, and then applying to 
pull back to the domain X and function /. By restricting the middle space to 
a finite graph, these constructions enjoy the simple topology of graph continua 
m- Monotonicity of p makes the theory oblivious to the many complexities of 
continua. 

We use the light factor A to measure length along the edges of M : For points a, b 
on an edge, the length between them is |A(a) — A(6)|. Topological total variation 
TTV (/) is the sum of all M’s edge lengths. 

A varilet basis is a finite collection {gi \ i = l...n} of real-valued piecewise 
monotone functions, called varilets, such that every linear combination ^ OiPi 
(each Oi G K) has topological total variation TTVi^Oigi) = ^ 

ilet basis is normal and independent in the sense that each TTV{gi) = 1 and 
I] a^gi = I] hgi only when each Oi =bi. 

The name “varilet” reflects the relationship of the basis functions to topological 
total variation. 

A varilet transform for piecewise monotone f is a varilet basis such that / is a 
linear combination: / = X) The Varilet Transform Algorithm is stated as 
a mathematical algorithm from which computational methods may be derived. 

Varilet filters for / are created by varying the coefficients, yielding filtered func¬ 
tions f' = Varilet filters manipulate topological total variation analo¬ 

gously to linear filters’ manipulation of energy. 


1.2 Related Work 

Varilets fit into the larger context of computational topology and data analysis 
mM- Although there may exist mathematical connections to persistent ho¬ 
mology [T^ and discrete Morse theory [T3], this paper focuses instead on the 
monotone-light factorization m, which for piecewise monotone functions may 
be seen as a decorated version of the Reeb graph m- Sometimes called the 
contour tree [S], and within persistence theory the merge tree m, the Reeb 
graph has been often used for simplification of scalar fields mm and has been 
exploited throughout computational topology. 

The present paper has connections to the work of Bauer et al. on persistence and 
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total variation [5] , and to work relating topological analysis to signal processing, 
including Guillemard et al. imiii] and Bauer et al. There are similarities 
in intent, but not in formalism, to Robinson’s notion of topological filter [22] . 

Our definition of topological total variation is similar to graph total variation as 
defined by Bresson [33] . Our definition agrees with the usual ID definition of to¬ 
tal variation, but conflicts with most multidimensional definitions |15] , including 
that used in image processing [23] , 

1.3 Contributions 

This paper introduces the varilet transform, an additive decomposition of scalar 
fields by independent normalized summands, which also additively decomposes 
a generalized total variation measure. 

The Varilet Transform Algorithm provides a mathematical skeleton for compu¬ 
tational methods. The algorithm is proved correct. 

Filtering / by varying the coefficients of its varilet transform provides an gen¬ 
eralization of simplification of scalar fields IHEl]. 

1.4 Organization 

This paper has two main sections: Section provides definitions, culminating 
in the Varilet Transform Algorithm. Section validates the algorithm, proving 
that it computes a varilet basis. 


2 Definition of the Varilet Transform Algorithm 

Following some preliminaries, we define piecewise monotonicity, topological total 
variation, and then varilet bases, transforms and filters. Finally, we state the 
Varilet Transform Algorithm. 

2.1 Preliminaries 

All spaces in this paper are compact metric spaces. A continuum is a con¬ 
nected component of a compact metric space; we assume all continua are non¬ 
degenerate (no isolated points). Continua include line segments, disks, spheres, 
simplexes, graphs having one-dimensional topology on their arcs (e.g. Reeb 
graphs), compact manifolds, as well the result of (appropriately) attaching to¬ 
gether other continua. 

We specify /’s monotone-light factorization, abbreviated “m.l.f.”, by simply 
listing fiMX, where is the monotone factor, M is the middle space, and A is 
the light factor. 
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For continuous / on compact metric space X, the middle space M has finitely 
many connected components; they are in 1-1 correspondence with X’s compo¬ 
nents. There are no relations among components of X, nor among components 
of M. 

f’s middle space M is a compact metric space; therefore assertions proven for / 
also apply to A. The m.l.f. of A is IMA, where 1 represents the identity function 
on M. 

We refer to the connected components of f~^y as contours. 

For any subset ^ of a topological space, we denote the interior by S°, the closure 
by S, the boundary by dS, and the complement by S'^. 


2.2 Piecewise Monotone Functions 


This section introduces piecewise monotone functions, which will be used through¬ 
out. They correspond in principle to use of tame functions [Mill]. 

The monotone-light factorization provides the basis for piecewise monotonicity. 
We give a general definition, followed by specialization to real-valued functions. 

Definition 2.1 (Piecewise Monotone Function). Suppose f : X ^ Y has m.l.f. 
fxMX. We say that A is locally monotone at p £ M whenp has a neighbourhood 
upon which A is monotone. Let M* denote the set of all points of M at which 
A is locally monotone. Then f is piecewise monotone when: 

(1) M* is dense in M; 

(2) M* has finitely many components; and 

(3) A is monotone on the closure of each component of M*. 


The closures of the components of M* are the monotone pieces referred to in 
the name, which we abbreviate as “p.m.”. 


The light factor A is a homeomorphism on each monotone piece. 


monotone. 


For p.m. scalar field / : X —)■ K, each monotone piece of the middle space M 
is a closed, non-degenerate real interval. Some interval end points are shared 
between two or more - but only finitely many - intervals; these are the points of 
definition |2.1| at which A fails to be locally monotone. Some interval endpoints 
are not shared - in this case A is locally monotone at the endpoint. The mono¬ 
tone pieces form a finite graph; the interval endpoints are the graph vertices; 
each interval defines a single edge. M may have more than one edge between a 
pair of vertices, and multiple edges may constitute a loop, but M does not have 
an edge connecting a vertex to itself. 

The middle space M of p.m. / : X —K is graph-theoretically and topologically 
identical to f’s Reeb graph [3T]. We will use both its graph structure and its 
point-set topology, but we will not be concerned with homotopy. When compact 
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metric space X has connected components Xi... then the middle space M 
comprises n disjoint graphs Mi... Mn, and each f\Xi is piecewise monotone 
with middle space Mi. 

It is well known that when X is simply connected then M is acyclic; when 
/ : [0 1] —M then M is linked chain of closed intervals. Graph continua appear 
in the literature [13 E!, but have not been previously utilized in relation to 
monotone-light factorization. 

For f : X ^ M., the light factor A is numerically strictly monotone along each 
edge of M. For an edge having endpoint at vertex V, we use A to characterize 
the edge as increasing or decreasing at V. If vertex V has both an increasing and 
a decreasing edge, then it is a saddle; in this case V terminates three or more 
edges. A vertex for which all edges have the same direction is an extremum, 
either a maximum or minimum. The extrema and saddles of M are called 
critical points. (This usage is more general than the classical notion, in the 
same spirit as homological critical values in [lOj.'l 

We use A to measure length along the edges of M: For any two points a,b on 
an edge (including endpoints), the length between them is |A(a) — A(6)|. 

When / is piecewise monotone, then so is its light factor A. 

2.3 Topological Total Variation 

Total variation TV (/) for differentiable / : [0 1] —>■ K is given by TV (/) = J \f'\. 
On multidimensional domain X C K”, total variation has a gradient formulation 
TV{f) = / |V/|. For non-differentiable functions, see definitions, examples and 
historical discussion in |15] . 

We provide an alternative definition: 

Definition 2.2. For p.m. / : V —>■ K having m.l.f. pMX, f’s topological total 
variation TTV (/) is the sum of all M’s edge lengths. 

TTV{f) = TV if) when / : [0 1] ^ K. 

However, when X is multidimensional then TTV and TV do not agree, seen 
as follows. The co-area formula m in equation Q expresses total variation as 
the integral of level set perimeter lengths Per{f,y), 

/ OO 

Per{f,y)dy. (1) 

-OO 

Whereas TTV is invariant under self-homeomorphisms of X, equation Q indi¬ 
cates that multidimensional TV is not; for example, consider a homeomorphism 
of X that stretches f’s level set perimeters. 

Topological total variation avoids dependence on the particulars of f’s domain 
X by measuring instead on f’s middle space M. TTV{X) = TTV{f). 
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When compact metric space X has connected components Xi...Xm then 
TTV{f) = E TTV{f\Xi). 


2.4 Varilet Basis 

Definition 2.3. A finite collection {gi :X— | z = l...n} of p.m. functions 
is a varilet basis when TTV ^i9i) — S l®*l choices of at G K. 

The functions pi are called varilets. A varilet basis is normal and independent 
in the sense that each TTV{gi) = 1, and TTV{J2o-i9i) = 0 iff each all Oi = 0, 
implying that every linear combination is unique. 


2.5 Varilet Transforms 

Definition 2.4. A varilet transform for p.m. f : X —>■ M. is a varilet basis 
{pi I i = 1... n} such that there exist positive amplitudes {at \ i = 1... n} 
giving f = J2 oii9i ■ 

Unlike the Fourier transform, there is no unique varilet transform for a piecewise 
monotone function; instead there are multiple varilet transforms, in this respect 
similar to wavelet transforms. 

Section [2.7| provides a mathematical algorithm having two input parameters: the 
function /, and a choice of lens C. The Varilet Transform Algorithm outputs a 
varilet basis {pi | z = 1.. .n} and amplitudes {ai | z = 1.. .n}. Sectionwill 
validate the algorithm by developing some elementary analytic topology. 

The Varilet Transform Algorithm provides one way to compute varilet trans¬ 
forms, but definition |2.4| may admit other methods as well. 

2.6 Varilet Filters 

Definition 2.5. Let {pi | z = 1... rz} be the varilet basis resulting from a varilet 
transform for f. Every choice of varilet filter coefficients {oi S M| z = 1... zz} 
defines a filtered version of f,f' = Yl aiPi. 

Since | z = 1... n} is a varilet basis, it follows that TTV (/') = ^ |ai|. 

2.7 Varilet Transform Algorithm 

This section introduces the mathematical objects and algorithm for computing 
varilet transforms. Objects introduced are the varilet lens C, and varilet supports 
E. The lens C plays the role of a parameter in the algorithm. 


2.7.1 Varilet lens 


The Varilet Transform Algorithm analyzes a function, the parameter /. This 
section introduces the only additional parameter, the transform’s lens C, so- 
called because it determines the selection and resolution of f’s analysis in rela¬ 
tion to critical points in /’s middle space M. 

We do not specify how to choose the lens C, but instead take an axiomatic ap¬ 
proach, relating specialized properties of C to resulting properties of the varilet 
transform. Looking forward to those results, the varilet transform’s twin addi¬ 
tive decomposition of function / and topological total variation TTV{f) is a 
universal property of the varilet transform, i.e. it holds for all C. Whereas, in 
the contexts of image segmentation, simplihcation and fractal analysis [1], we 
will require lens C to have special properties. 

Definition 2.6 (Varilet Lens). A subset C C M is a constant-boundary region 
when C is closed, connected, has nonempty interior, and A is constant on dC. 

Let I be a finite sequence of natural or real numbers that will serve as indices. 
An indexed collection C = {Ci | i G 1} o/ constant-boundary regions is nested 
when indices i < j imply either Ci 3 Cj or CiCiCj =0. 

A varilet lens for / is nonempty nested collection C of constant-boundary regions 
such that IJCi = M. 

I 

The root regions of varilet lens C are the connected components of M, and the 
root indices are the indices of the root regions. 

For each i G I, we define Ci ’s successors Si be the collection of all maximal 
(by inclusion) Cj C Ci. When Cj is a successor of Ci, we say that Ci is the 
predecessor ofCj. 

When unambiguous, we suppress the index set I, writing {Ci ...} instead. 

The successor relation exactly reflects the tree of inclusions within C. 

The root regions are the only constant-boundary regions having empty bound¬ 
ary. The root indices comprise an initial subsequence of T. 

The nesting structure of C suggests a multiresolution lens, with fine resolution 
for tightly nested sets in C and coarse resolution for sparsely nested sets. 

2.7.2 Varilet Supports 

From the varilet lens C we derive the varilet supports F, which play a central 
role in the Varilet Transform Algorithm. 

Definition 2.7. Let C = {Ci \ i (il) be a varilet lens for f. Then C’s varilet 
supports comprise the collection V = {Di \ i G I{, where each varilet support 
is defined as Di = {Ci \ [_)Cj). 

j>i 
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The name “support” is discussed in section [2.7.3 
Transform Algorithm. 


after introducing the Varilet 


Due to the simplicity of the finite-graph continuum topology of /’s middle space 
M, we can state some elementary properties: 

• A varilet support Di is not necessarily connected, but has only finitely 

many components, each of which is the closure of a nonempty connected 
interior. 

• Each component of each Di is a connected graph fragment cut out of M; 

these fragments overlap only at the cut points, which may be either regular 
or critical points. 

• D covers M, and supports can intersect only at their boundaries. 

• Each Ci = U{Dj | Dj C Ci} and each dCi C U{dDj \ Dj C Ci}. 

• When no Cj C Ci, then Di = Ci- 


2.7.3 Varilet Transform Algorithm 

The Varilet Transform Algorithm is a mathematical algorithm resulting in a 
varilet transform for / based on lens C. Since there are many possible C, the 
algorithm can produce as many different varilet transforms for /. 

The output of the Varilet Transform Algorithm is an indexed collection of func¬ 
tions {gi : X ^ R \ i G I}, and a correspondingly indexed collection of positive 
reals {ai \ i € I}. 

We equate the collection of functions to the varilet transform of f with lens C, 
written: 

= hei}, 

and we call {ai | i G 1} the transform’s amplitudes. 

Definition 2.8 (Varilet Transform Algorithm). Suppose p.m. f : X —>■ M. has 
m.l.f pMX. Let C = {Ci | i S 1} be a varilet lens for f, and let D = {Di ...} 
be the varilet supports. 

Each function gi will be created by first creating functions Xi,"fi : M —> K, then 
defining gi on X via f’s monotone factor: gi = jiO p,. 

For each i G I, in any order, or asynchronously, define gi and ai: 


(1) Define ai = TTV{X\Di). 


Note: Lemma 3.4 will show X\Di piecewise monotone. 


(2) Let M* be the component of M containing Di. 


Define Xi on M* as the unique extension of X\Di that is constant on 
each component of M* \ Di. 

Define Xi = 0 on all other components of M. 
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Note: Lemma 3.6 will show that there exists a unique piecewise constant 


extension. Lemma 3.4 will show piecewise monotone and TTV{Xi) = 


(3) Define 7 ^ on M* hy scaling and shifting A^. 
VpGM*, 


, -1 


7i(p) = 


A*(p) 


when i is a root index 


Oj {Xi{p) — X{dCi)) otherwise. 


Define ji = 0 on all other components of M. 


Note: Lemma 3.4 will show 7 ^ piecewise monotone and TTV ( 7 ^ 


(4) Define gi = o p. 


Note: Lemma 3.4 will show gi piecewise monotone and TTV{gi) = 1. 


The Varilet Transform Algorithm produces functions for which we have not yet 
proved any properties; in fact, we must prove that the algorithm creates well- 
defined functions. The next section proves that y{f,C) = {gi | i S 1} is a 
varilet basis expressing / = X) 0 ^* 5 * ■ 

We can now better motivate the name “support”. Each Di supports functions 
Xi (and 7 i) in the somewhat nonstandard sense that Xi is nowhere constant on 
Di but is constant on each component oi M \ Di. “Support” usually means 
“non-zero” but here means “non-constant”. 
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3 Validation of the Varilet Transform Algorithm 


In this section we prove: 

Theorem 3.1 (Varilet Transform Theorem). The result of the Varilet Trans¬ 
form Algorithm, 'V{f,C) = {gi \ i € X}, is a varilet basis, and f = 


The rather involved proof compensates for a lack of preexisting mathematics 
by creating a small theory having its own vocabulary. Almost all work uses the 
middle space and light factor as surrogates for domain and function, using the 
monotone factor to connect the two when required. 


Section 

Section 


3.2 


3.3 


establishes that the Varilet Transform Algorithm is well-defined, 
shows that the algorithm additively decomposes / and TTV{f). 


Section 3.4 shows that the algorithm creates is a varilet basis, completing proof 
of theorem 13.11 


Throughout, we work with p.m. / : V —>■ K having m.l.f. fiMX. 


3.1 Flat Extensions and Constant-Boundary Functions 

We start by identifying two classes of functions on /’s middle space M. 

Definition 3.2 (Flat Extension). Consider any closed D C M and any function 
TT : — )■ K, and letir* : M — > K be an extension o/tt. Then tt* is a flat extension 
when TT* is constant on each component of M \ D. 

Varilet lens C causes the varilet transform and filtered versions of / to be similar, 
in the sense that they all have form tt o p,, where tt is the following type of 
function. 

Definition 3.3 (Constant-Boundary Functions), tt : M ^ M. is a constant¬ 
boundary function for C when tt is constant on each dCi, for j G X. 

A is constant-boundary for C. 

In this paper we will define new constant-boundary functions tt, identifying 
useful relationships between functions / and f' = tt o p, where tt is substituted 
for A in f’s monotone-light factorization. 

C’s constant-boundary functions are closed under composition with any contin¬ 
uous function C : K” —>■ R, i-e. Tr(p) = ■ ■ ■'^n{p)) is constant-boundary, 

for any choices of constant-boundary tti ... 7r„. 

In the Varilet Transform Algorithm, once we have established that the functions 
Ai, 7 i are well-defined, it will follow that each is constant-boundary for C. 


3.2 Varilet Transform Algorithm Well-Defined 

In this section we validate the notes attached to the Varilet Transform Algo¬ 
rithm, thereby establishing that each function Ai, 7 ^ and gi is well-defined, p.m.. 
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and has topological total variation as indicated in the statement of the algo¬ 
rithm. 

Section [3.2.1| proves elementary properties of restrictions and extensions to piece- 
wise monotone functions. Section |3.2.2| proves existence and uniqueness of flat 
extensions to restrictions of constant-boundary functions, completing the proof 
that the Varilet Transform Algorithm is well-defined. 


3.2.1 Restrictions & Extensions of Piecewise Monotone Functions 


Statement (a) of the following Restriction-Extension Lemma validates step (1) 
the Varilet Transform Algorithm. Once lemma 3.6 establishes that the functions 
Ai, 7 i and gt are well defined, it will follow from Restriction-Extension Lemma 
(b) & (c) that each is p.m. and has topological total variation as indicated in 
the algorithm, validating steps (3) & (4) of the algorithm. 


Lemma 3.4 (Restriction-Extension). Choose any closed D <Z M having finitely 
many components each having nonempty interior, and let the restriction f = 

(a) f is p.m. and has m.l.f. p!M'\' where p,' = p\{p ^D), M' = D, and 
X' = X\D. 

(b) TTV (/') is equal to the sum of edge lengths in D, noting that edge lengths 
in M and M' are determined by X and X respectively. 

(c) Suppose g is a flat extension of /'; then g is piecewise monotone and 
TTV{g) = TTVif). 


Proof. Because p is monotone, p~^D has components in 1-1 correspondence 
with Zl’s components and is therefore a compact subspace, and thus /' has 
a monotone-light factorization. Every contour of f' is a contour of /, so it 
follows that M' is homeomorphic to D, and so p' = p\{p~^D) and A' = X\D 
by uniqueness of monotone-light factorization, f' is p.m. since M' is a finite 
graph. When considered as the middle space of /', D has vertices at what were 
cut points of edges of M. These observations prove (a) & (b). 

We now show (c). Let pgMgXg be g's m.l.f. For every component K oi M \ D, 
g is constant on p~^K. Therefore pg maps p~^K to a point in Mg. Therefore, 
Mg consists of the components of D°, with some combination of boundary 
points glued together. This results in a finite graph and therefore g is piecewise 
monotone. Mg is a quotient of M; the quotient map is a homeomorphism on 
D°, from which it follows that Xg measures edge length identically to A, and 
thus TTV{g) = TTV{f'). □ 


The following Lemma provides a recipe for computing topological total varia¬ 
tion. 

Lemma 3.5 [TTV Decomposition Lemma). Let Di, i = 1.. .n, he any collec¬ 
tion of compact subsets covering M that may intersect only at their boundaries. 
Then TTV[f) = ^ TTV[X\Di). 
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The TTV Decomposition Lemma implies TTV{f) = where ai are the 

amplitudes output by the Varilet Transform Algorithm. We will use the lemma 
again in the proof of the Varilet Transform Theorem, in section [3. 4. 5[ 


3.2.2 Flat Extension Lemma 


The following lemma validates step (2) of the Varilet Transform Algorithm, 
thereby completing the proof that the Varilet Transform Algorithm is well- 
defined. 

Lemma 3.6 (Flat Extension Lemma). Choose a varilet lens C = {Ci \ i G 1} 
of /, and let V = {Di ...} be the varilet supports. Choose any tt : M —>■ K from 
C’s constant-boundary functions. 

Then for every varilet support Di, i & dd: 

(a) There exists a flat extension tt* : M —>■ K of TT\Di. 

(b) When M is connected then it* is unique. 

(c) TT* is a constant-boundary function for C. 


We will need: 

Lemma 3.7. With C, D and tt as in the previous lemma: 

(a) TT is constant on the boundary of each component of M \ Di. 

(b) Let Zi ... Zn he the unique values among Tr{dCi) and all TT{dCj), for suc¬ 
cessors Cj € Si. Then the only values that tt may have on then boundary 
of any component of M \ Di is one of Zi... Zn, and for each zj- there exists 
at least one component having boundary upon which tt has this value. 


Proof Let K = M \ D, = CfU (UCj). Then dK = dC, U {CdCfl. 

Si Si 

Ci’s successors Cj G Si are disjoint, but a successor Cj may have boundary 
intersecting dCi. Because each is connected, each successor Cj lies in a single 
component of K. To prove (a), we must show that when Ci’s successors Cj, Ck 
lie in the same component of K, then TT(dC'j) = ir^dCk). But when they lie 
in the same component, then there must exist a path connecting them in Cf, 
from which it follows that TT^dCi) = TT(dCj) = Tr{dCk). Statement (b) follows 
immediately. □ 


Proof of Flat Extension Lemma (3.6). Using the notation of the preceding proof, 
extend 7r|Di to tt* : M —>■ M by defining tt* on each component ot K' C K to 
have constant value TT{dK'), thereby proving (a). Statements (b) & (c) follow 
immediately. □ 
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3.3 Additive Decomposition 

In this section we prove part of the Varilet Transform Theorem (3.1). 
Lemma 3.8 (Additive Decomposition). / = 


To prove the Additive Decomposition Lemma, we first link constant-boundary 
function values in a tree structure, and then use an inductive argument exploit¬ 
ing the links. The Link Lemma is proved in the next subsection, with proof of 
the Additive Decomposition Lemma in the section following. 


3.3.1 Link Lemma 

The varilet supports T) are fully determined by choice of varilet lens C. The 
inclusions within C impart an identical tree structure on V. The next lemma 
describes how this tree is manifested in the boundaries of the varilet supports 
in (D, and its implications for constant-boundary functions. 

Lemma 3.9 (Link Lemma). Choose a varilet lens C = {Ci \ i € X}, and let 
V = {Di ...} be the varilet supports. Choose any index i g2. 

Then for each successor Cj € Si there exist points pi G dDi, qj G dDj, such 
that for any constant-boundary function it for C, '^{pi) = T^iqj) = TT{dCj). 

We call {pi,qj) a link pair, because it associates a point in predecessor Di to 
a point in successor Dj having equal 7r-value, independent of the choice of 
constant-boundary tt. There exists a link pair (j)i,qj) for all indices j G I 
except the root indices. 

Proof. Recalling definition |2. 7. 2[ it follows from pairwise disjointness of the Cfs 
successors that dCj intersects both dDi and dDj. □ 


3.3.2 Proof of Additive Decomposition 


Using functions Xi,^i defined in steps (2) & (3) of the Varilet Transform Algo¬ 
rithm, we will use the Link Lemma to show A = ^ OiTi- Then, using p, to pull 
back from M to /’s domain X, we get f = C(i9i- 

We use the following: 

Lemma 3.10 (Zero Varilet Lemma). Choose any varilet support Di. 

Then jilDj = 0 for all indices j such that Cj (f Ci. 


Proof. Step (3) of the Varilet Transform Algorithm states that ^i{dCi) 


When Cj <fL Cj 


then Dj C Cf 


Since 7 ^ is constant on each such Dj, the 


Lemma (3.9) implies ')i{dDj) = ')i{dCi) = 0. 


= 0 . 

Link 

□ 


We now continue with: 
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Proof of Additive Decomposition Lemma (3.8). For each index z S I let 


TT, = and Ei = [jDj. (2) 

j<i j<i 

We will show that TTi\Ei = X\Ei. Letting n = max(I), since 7r„ = A and 
En = M, this will prove A = ^ oti'yi- 

We proceed by induction on increasing indices i G X. 

When z is a root index, steps (2) & (3) of the Varilet Transform Algorithm state 
that 'Ki\Di = ai"fi\Di = X\Di. All indices j < z are also root, so TTi\Ei = X\Ei. 

When z is not a root index, we can assume by induction that T^j\Ej = X\Ej for 
all indices j < i. We must show 


X\Ei = ni\Ei. 


By the Zero Varilet Lemma (3.10), jilDj = 0 for all j < i. Let k be the index of 
Cfs predecessor. For index j, k < j < i, Cj and Ci are disjoint, so = 0, 


again by the Zero Varilet Lemma (3.10). Therefore, for any j, k < j < i, 


n^\Ej = TTklEj 

= X\Ej by induction. 


Thus, it suffices to show 


X\Di = ni\Di. 


Expanding the definition of tt^. 


j<i 

= ai-fi\Di + TTfclA, 


and thus it sufhces with 


A|A = aijt\Di +7rfe|A- 


(3) 


TTfe is constant on Di] we can identify the constant value: The Link Lemma 
provides linked pair {pk,qi), where pk G dDk, qi G dDi, and nkipk) = T^k(.qi) = 
TTkidCi). 

But by induction, Trk{dCi) = X{dCi). Thus equation (§ is equivalent to 

A|A = aailDi + X{dCi). 


The definition of 7^ in step (3) of the Varilet Transform Algorithm completes 
the proof. □ 
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3.4 Varilet Basis 


Proof of the Varilet Transform Theorem is completed in this section. 
Lemma 3.11 (Varilet Basis). y{f,C) = {gi \ i G 1} is a varilet basis. 


We start by defining link recursion, a definitional schema whereby the Link 
Lemma (3.9) enables piecewise definition of constant-boundary functions for C. 
Section |3.4.2| uses link recursion to define the filter factor ip, proving that it 
describes varilet filtered versions of /. Section [3.4.3] describes how proof of the 
Valrlet Basis Lemma (3.11) can be split into two cases A & B, and we prove 
case A. The development continues in section 3.4.4 with dehnition of the filter 
quotient f, and identification of the m.l.f. of varilet filtered functions, thereby 
enabling proof of case B, in section [3. 4. 5[ 


3.4.1 Link Recursion 


This section describes link recursion, a method for piecewise definition of tt : 
M —>■ M, a new constant-boundary function for C. Link recursion is based on 
the Link Lemma (3.9). 


To define a constant-boundary function tt : M — >■ M using link recursion, we 
iterate through the indices i G I in increasing order, defining tt on Di at each 
iteration. 


The root indices i come first in the iteration; for each root index i, TTijUj may 
be independently defined as any p.m. function that is constant on each dCj, for 
Cfs successors Cj. 

Following the root indices, for each iteration i, let k < i be the index of Cfs 
predecessor, and let {pk,gi) be their link pair. Tr\Di may be independently 
defined as any p.m. function such that tt is constant on each dCj, for Cfs 
successors Cj, and such that tt is constant on dCi with value TT{dCi) = TT{pk), 
noting that the value of T^fpk) bas previously been defined in iteration k. 

Link recursion results in a well-defined, continuous, p.m. function tt, a constant¬ 
boundary function for C. 


3.4.2 Varilet Filter Factor 

In this section we discuss varilet filters in relation to the Varilet Transform 
Algorithm. 

The Varilet Transform Algorithm results in functions y{f,C) = {gi...} and 
amplitudes {at...} such that f = For filter coefficients A = {ai\i G I}, 

we define notation for the varilet filtered function: 

= ^a^gi. 
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To prove that is a varilet basis, we need to show that TTV{J2°-i9i) — 

^ \ai\, for every choice of {ui \ i G I}. Therefore it should not be surprising 
that proof of the Varilet Basis Lemma (3.11) uses results that apply to varilet 
filters generally. 


Definition 3.12 (Varilet Filter Factor ■0). Choose a varilet lens C = {Ci \ i G 
T}, and let T) = {Di ...} be the varilet supports. 

For any choice of filter coefficients A = ...}, the filter factor : M —>■ K 

is defined by link recursion: 


'fpAip) 


{ui/ai)X{p) p G Di, when i is a root index; 

{oi/ai){\[p) — X{dCi)) + p G Di, iterating over increasing i. 


is a p.m. constant-boundary function for C. 

Let F be the collection of all indices j such that aj = 0; then tpji is a flat 
extension oi D = iPaK U Dj). 

The name “filter factor” is motivated by the substitution of ip for A in f’s 
monotone-light factorization in: 

Lemma 3.13 (Filter Factor). Let y{f,C) = {gi \ i & D}, let A= {oi...} he 
any choice of varilet filter coeff dents, and let ip a be the filter factor. 

Then OiPi = ipAO p. 


Proof. Let ip = ip a, let /' = ip o p, and let Xi = p ^Di for each i G T. We 
proceed by induction on increasing indices i gX. 


Suppose z G I is a root index; the definition (3.12) for ip together with root- 
index definition of 7^ in step (3) of the Varilet Transform Algorithm imply that 
f'\Xi = aigi\Xi. Any j < i are also root indices; therefore /' = 

.ux,. 

3<'b 

Now suppose z € Z is not a root index. We may inductively assume 

/' = on LJ^J- 

j<i j<i 

We must show 

/' = Yl°-j93 on [jXj. 

j<i j<i 

By the Zero Varilet Lemma ( 3.10[ ), gi\Xj = 0 for all j < i, and so by induction 
it suffices to show 

/' = Ylaj93 on Xi. (4) 

j<i 

Let k < i he the index of Cfis predecessor. gj\Xi is constant for each j < k, 
and therefore so is their sum. Using the Link Lemma (3.9) we can identify the 
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constant value: 


Y,a,ig,\X,) = f'idX,) 

j<k 

= i’idX,) 


by induction. 


For index j, k < j < i, Cj and Ci are disjoint, and so gj\Xi = 0 by the 
Zero Varilet Lemma (3.101. Therefore, expanding the sum in equation ([4]), it is 
sufficient to show 

f = arQi + i^idXi) on X,, 


which follows from definition of "0 and the non-root-index definition of 7 ^ in step 
(3) of the Varilet Transform Algorithm. □ 


3.4.3 Proof of Varilet Transform Theorem, Part A 


We must show that X{f,C) = {gi...} is a varilet basis; i.e. for every {a^...} 

TrF(^a,5,)=El“*l- 

In light of the Filter Factor Lemma ( |3.13[ ), this is equivalent to 

TTVi^og)=J2\a.\. (5) 


The intuition for the proof is that the filter factor '0 multiplicatively stretches 
the lengths of the edges in each varilet support Di, using coefficient as stretch 
factor. To measure topological total variation oi f = ip o g = ^ aigi we must 
identify /'’s middle space and light factor. When the coefficients at are all 
nonzero, then the m.l.f. of f is g,Mip, from which we easily show ©• A small 
complexity enters the proof when one or more are zero; in this case the middle 
space of /' is not M. We identify the middle space as a certain quotient of M; 
this will allow us to show ([^. 


We break up the proof of the Varilet Transform Theorem (3.1 1 into two parts: 
Part A covers the case where all filter coefficients at 7 ^ 0; part B covers the case 
where one or more = 0 . 


Proof of theorem 3.1 , part A (oi ^ 0 for all i Gl). Let g!M'X' denote the m.l.f. 
of f = = 0 o Filter factor ip is light because A is and each 7 ^ 0, 

and therefore by uniqueness of monotone-light factorization g' = g, M' = M, 
and A' = ip. 

Because M' = M, C is a varilet lens for /'. The TTV Decomposition Lemma 
(3.5) states that TTV{f') = TTV{ip\Di). By definition (3.12) of ip, 

Y, TTViiPm = ^(|a,|M) FTF(A|A) 

= Y 

completing the proof. □ 
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3.4.4 Varilet Filter Quotient 


The section identifies the monotone-light factorization of filtered function /' = 


Lemma 3.14 (Filter Quotient Lemma). Suppose p.m. f has m.l.f. p,MX. Choose 
a lens C = {Ci \ i € I}, and choose filter coefficients A = {oi.. .}, and suppose 
varilet filtered f = ^{f,C,A) has m.l.f. p!M'\'. 

Then: 

(a) M' is a certain quotient of M, described in the proof. 

(b) Denote the quotient map as (j)j\ : M —>■ M'. Then is monotone. 

(c) jj! = (/)^ o pL, M' = <f>j^{M), and A' = o for any right inverse (jf^. 


Definition 3.15 (Filter Quotient f). 
called Al’s filter quotient. 


The quotient map in lemma 3. If 


IS 


Proof of lemma 3. If The statements pertaining to the case where all Oi ^ 0, 


and therefore M' = M, follow from the arguments in the proof of part A of 
Varilet Transform Theorem in the previous section. We prove the lemma for 
the case that one or more = 0. Let filter factor ip = "0^. 


Let P be the nonempty collection of all indices j € D for which Oj = 0, and 
assume I \ is nonempty (since otherwise ^ Oigt = 0). 

We claim that M' is the quotient of M that identifies all points in each compo¬ 
nent of ( U Dj). Let p = 0^ be the quotient map. 


We proceed by showing that p and any right inverse 
diagram commute, recalling that Filter Factor Lemma 
atgi = ipo p. 


0’" make the following 
(3.13) states that /' = 


X M M 



It is obvious that p{M) is a finite graph. We show that p{M) is the middle 
space M' of f by showing that 0 o p, is monotone and that 0 o 0’’ is light, for 
any choice of right inverse 0’'. 

po pL is monotone because p is monotone: For any point q G p{M) the set p~^q 
is connected, being either singleton or a connected component of ( U DA. 

We show that any right inverse 0’’ can be chosen. Let p be any point chosen 
from p~^q; we show that the value of p{p) does not depend on the choice. When 
p~^q is singleton then the choice is irrelevant. Otherwise, p~^q is a connected 
component of ( U Dj), and 0 is constant on each such component. 
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To see that ij) o (jf is light, choose any q G and let z = {ip o (j)''){q). We 

show that (V' o (jf)~^z is a finite collection of points. Because (^{M) is a finite 
graph, if {ilio(jf)~^z were infinite then it would contain an open set of 4>{M), in 
which case 'tp~^z must contain an open set of M. But any open set of M upon 
which ip is constant must lie in ( U Dj)^ which is mapped by (/) to a finite set 

of points in <p{M), a contradiction. 

We have now confirmed that diagram § commutes, i.e. that the m.l.f. of f' is 
fi' = (po fx, M' = (p{M), and X' = ip o cp''. □ 

We state some immediate consequences regarding the filter quotient <p'. 
Lemma 3.16 (Filter Quotient Lemma). 

(a) When Ui ^ 0 then: 

- (p is a homeomorphism on D°. 

- (p{Di) may have fewer components than did Di, but may not have more. 

- (p~^{d(p{Di)) C dDi, but (p may also map points in dDi to (p{Di)°. 

(b) The collection {(p[Di) | ^ 0} covers M' and the sets (p{Di) intersect 

only at their boundaries. 


3.4.5 Proof of Varilet Transform Theorem, Part B 


It is now easy to complete the proof of the Varilet Transform Theorem (3.1). 


Proof of theorem 3.T part B. Let /' = .^(/, C, A)] we must show that TTV (/') = 

EkI- 

The Filter Quotient Lemma (3.16) and the TTV Decomposition Lemma (3.5) 
allow us to compute TTV{f'), 


TTVif) = Y, TTViX'lcPiD,)) 


(7) 


When Oi 0, since cp is a homeomorphism on D°, it follows from the Filter 
Quotient Lemma (3.14) and definition (3.12) of filter factor ip that 


TTV{X'\PiD,)) = TTViiPm 

= (|a,|/a,) TrV(A|A) 
= |a^|. 


Thus, equation ([^ becomes 


TTVif) 


ai^O 




completing the proof. 


□ 
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